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Abstract 



The paper presents a study of a solar dynamo model operating in the bulk 
of the convection zone with the toroidal magnetic field fiux concentrated in the 
subsurface rotational shear layer. We explore how this type of dynamo may 
depend on spatial variations of turbulent parameters and on the differential ro- 
tation near the surface. The mean-field dynamo model takes into account the 
evolution of magnetic helicity and describes its nonlinear feedback on the gen- 
eration of large-scale magnetic field by the a-effect. We compare the magnetic 
cycle characteristics predicted by the model, including the cycle asymmetry (as- 
sociated with the growth and decay times) and the duration - amplitude relation 
(Waldmeier's effects), with the observed sunspot cycle properties. We show that 
the model qualitatively reproduces the basic properties of the solar cycles. 



1. Introduction 

The sunspot 's activity is organized on large scales, forming the Maunder butterfiy 
diagram. It is believed to represent the time - latitude pattern of the large-scale toroidal 
magnetic field generated in the convection zone. Another component of the solar activity is 
represented by the global poloidal magnetic field extending outside the Sun and shaping the 
solar corona. Both components synchronously evolve as the solar 11-year cycle progresses. 
The global poloidal field reverses the sign in the polar regions near the time of maxima 
of sunspot activity. Most of the current solar dynamo models suggest that the toroidal 
magnetic field that emerges on the surface and forms sunspots is generated near the bottom 
of the convection zone, in the tachocline or just beneath it in a convection overshoot layer 



(see, e.g.. 



Ruediger & Brandenburg 


1995; 


Choudhuri et al. 


1995 i 


Tobias & Weiss 


2007 
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Parker|[l993l ). The belief in a deep-seated solar dynamo comes from the fact that this region 
is sufficiently stable to store magnetic flux despite the magnetic flux-tube buoyancy effect. 
However, observations of rotation rates of emerging magnetic flux within the latitude bands 



±30° support a concept of relatively shallow sunspots ( Benevolenskaya et al.||1999 ) possibly 



rooted within the subsurface rotational shear layer. This concept has support from local 



helioseismology as well (Birch 2011) . 



There are further theoretical arguments that the subsurface angular velocity shear 
can play an important role in the dynamo process distributed in the convection zone 
( Brandenburg|[2005 ) . In our previous paper ( Pipin fc Kosovichev||2011 ) (PKll) we proposed 



a model of a subsurface-shear shaped solar dynamo. Our model shows that allowing the 
large-scale toroidal magnetic fleld to penetrate into the surface layers of the Sun changes the 
direction of the latitudinal migration of the toroidal fleld belts and produces the magnetic 
butterfly diagram in a good qualitative agreement with the solar cycle observations. The 
dynamo wave penetrates close to the surface and propagates along iso-surfaces of the angular 
velocity in the subsurface rotational shear layer in agreement with the Parker- Yoshimura 
rule ( Yoshimura||1975 ) . The standard boundary condition typically used in dynamo theories 



is to match the internal solution to the potential magnetic fleld extending outside of the 
dynamo region. This boundary condition does not allow to the toroidal component to 
penetrate to the surface. 

In our previous model the penetration of toroidal magnetic flelds to the surface was 
modeled by a special boundary condition at the top of the dynamo region. This boundary 
condition was formulated as a linear combination of vacuum (potential fleld) and perfectly 
conducting plasma conditions. The perfectly conducting part results in an increase of 
the toroidal component of the large-scale magnetic fleld at the boundary. Such boundary 
condition, used in PKll, models a partial penetration of the toroidal fleld into the solar 
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atmosphere, but from the physical point of view such formulation is rather artificial. The 
penetration to the surface can be modeled more physically by extending the computational 
domain close to the surface and using the magnetic diffusivity profile that follows from the 
standard solar interior model. This diffusivity decreases toward to the surface and results 
in increasing of the toroidal magnetic field in that direction (and an increase of the gradient 
of the toroidal magnetic field as well). The decrease of the turbulent diffusivity, tjt ^ -u'l^ 
(where u' is the convective RMS velocity and I is the mixing length) is predicted by the 
mixing-length theory of the solar interior. This motivates us to extend the integration 
domain from (0.71 -i- O.97)i?0, used in PKll, to (0.71 -i- O.99)i?0. The convection model 



of Stix (2002) predicts that towards to the surface the mixing-length, ^, decreases much 
faster than u' increases. Figure Ic shows the radial profile of the turbulent diffusivity in 
the convection zone model. In this paper we study how the sharp decrease of the magnetic 
diffusivity infiuences the strength and distribution of the toroidal field in the upper layers 
of the convection zone. 

There is another reason for extending the integration domain closer to the surface. In 
the near surface layers the density stratification gradient is very strong compared to the bulk 
of the convection zone. The mean-field theory predicts a downward turbulent drift of the 



large-scale magnetic filed in the presence of the density stratification gradient (Kitchatinov 



1991) (similar results were obtained by Pipin 2008). The effective downward drift of 



large-scale magnetic field results from magnetic fiuctuations in the stratified turbulence. 



It can be interpreted as follows (see Kitchatinov] 1991). The intensity of the magnetic 



fiuctuations 6^ = fiopu^'^ (p is the mean density) rises in the direction of the density gradient 

because the turbulent RMS velocity varies slower than p. Random Lorenz forces, which 

are induced by small-scale magnetic fields b and large-scale field B, produce fiuctuating 
(V X b) X B 



fiows ^ Tc. The resulted electromotive force x b is perpendicular to the 

MP 

large-scale field. This can be interpreted as an effective downward velocity drift of the 
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large-scale magnetic field ( Kitchatinov 1991). The theory also predicts that this kind of 



turbulent pumping is quenched by the infiuence of the Coriolis force, which results in the 
velocity of the effective drift to be greatest near the surface where the density stratification 
is strong. Thus, qualitatively, this effect works similarly to so-called "topological pumping" 



(Drobyshevski & Yuferev 1974) . 



Our study includes an equation of the magnetic helicity evolution, proposed by Kleeorin 



& Ruzmaikin (1982) and Kleeorin & Rogachevskii (1999). This equation describes the 



balance between the small-scale turbulent magnetic helicity and the large-scale magnetic 
helicity generated by the dynamo process, and has been used in many previous dynamo 



studies (e.g., Brandenburg & Subramanian (2005) and references therein). One of our 
goals is to explore nonlinear feedback of the magnetic helicity on the basic properties 
of sunspot cycles, e.g., the relationship between the rise and decay times, and between 
the length and strength of the cycles. The results for a dynamo model in a single-mode 



approximation (Kitiashvili & Kosovichev 2009, 2010) have suggested the importance of the 



nonlinear magnetic helicity effects for the solar-cycle behavior. The next section describes 
the formulation of the 2D mean-field dynamo model, including the basic assumptions, the 
reference model of the solar convection zone, and input parameters of the large-scale fiows. 
Section 3 presents the results and discussion. The main findings are discussed in Section 4. 



2. Basic equations 

The dynamo model is based on the standard mean-field induction equation in perfectly 
conductive media (Krause and Radler, 1980): 

— = V X (S+V X B) 

ot 
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where = u x b is the mean electromotive force, with u, b being the turbulent fluctuating 
velocity and magnetic fleld respectively; U is the mean velocity. General expression for E 



was obtained by Pipin (2008) (hereafter P08). Following Krause and Radler(1980) we write 
the expression for the mean electromotive force as follows: 



Ei = {aij + jij) Bj - rjijkVjBk. 



(1) 



Tensor aij represents the alpha effect, including the hydrodynamic and magnetic helicity 
contributions, aij = Ca^paW) ^i^^ Oa^f' + a[f^\ where the hydrodynamical part of the 



a-effect, aif\ and the quenching function, ip^^ are given in Appendix (see also in Pipin 



& Kosovichev 2011), the parameter Ca controls the amplitude of the a-effect. The 
hydrodynamic a-effect term is multiplied by sin^ 9 [6 is colatitude) to prevent the turbulent 
generation of magnetic fleld at the poles. The contribution of the small-scale magnetic 
helicity x = a-b (a is a fluctuating vector-potential of magnetic fleld) to the a-effect is 
deflned as a^^^ = C[f^x •> where coeflicient C\^^ depends on the turbulent properties and 
rotation, and is given in Appendix. The other parts of Eq.Q represent the effects of 
turbulent pumping, 7^^, and turbulent diffusion, r]ijk> They are the same as in PKll. We 
describe them in Appendix. 

The nonlinear feedback of the large-scale magnetic fleld to the a-effect is described as a 
combination of an "algebraic" quenching by function ip^. W) (see Appendix and PKll), and 
a dynamical quenching due to the magnetic helicity conservation constraint. The magnetic 



helicity, x •> subject to a conservation law, is described by the following anzatz (Kleeorin & 



Ruzmaikin 1982; Kleeorin & Rogachevskii 1999): 



(2) 



where Tc is a typical convection turnover time. Parameter controls the helicity 
dissipation rate without specifying the nature of the loss. Generally, we can expect that 
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the formulation of the the hehcity loss term in Eq.([2]) affects properties of the dynamo 



solutions. This is suggested by results that can be found in the literature (Brandenburg 



et al.||2007[ |Mitra et al.||2010[ [Guerrero et aL]|2010t [Mitra et al.||201lD . The physics of helicity 
loss is poorly understood, and influence the various processes of the helicity flux loss on 
the properties of magnetic cycles deserves a separate study. To reduce the number of free 
parameters in the model we consider the simplest form of helicity flux loss. The parameter 
controls the amount of the magnetic flux generated by the dynamo. This amount can be 
roughly estimated from observations. We use the range of that gives the total magnetic 
flux of the order of ^ 10^^ - 10^^ Mx in agreement with observations rtSchrijver & Harvey 



1994). Another parameter controlling the helicity dissipation in our model is Tc. It is given 



by the solar interior model. It seems to be reasonable that the helicity dissipation is most 
efl&cient in the near surface layers because of the strong decrease of Tc (see Figure lb). 



We use the solar convection zone model computed by Stix (2002), in which the 



mixing-length is deflned as ^ = aMLT |A(p)| , where A^*') = Vlog p is the pressure variation 



scale, and aMLT = 2. The turbulent diffusivity is parametrized in the form, tjt = Cj^r]). 



(0) 



where ry^^ = — is the characteristic mixing-length turbulent diffusivity, £ and are the 
typical correlation length and RMS convective velocity of turbulent flows, respectively and 
Cj^ is a constant to control the intensity of turbulent mixing. In the paper we use = 0.05. 
The differential rotation proflle, Q = Qofn (^^m), ^ = '^/Rq-> M = cos0 is a modifled version 



of an analytical approximation to helioseismology data, proposed by Antia et al. (1998), see 
Figure la. 

We use the standard boundary conditions to match the potential fleld outside and the 
perfect conductivity at the bottom boundary. As discussed above, the penetration of the 
toroidal magnetic fleld in to the near surface layers is controlled by the turbulent diffusivity 
and pumping effect (see Figures Ic and Id). 
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Model 








min 


G 


TDPolar 

G 


Aw 


Period, 
Yr 


PI 








-160 


1000 


7.2 


0.74 


13.4 


P2 


0.03 


- 


0.245 


-30 


450 


3.5 


0.94 


11.4 


P3 











60 


0.15 


0.89 


8.6 


Dl 

D2 


0.03 


- 


0.645 
0.091 


-160 


450 
1600 


3.6 

UA 


0.78 
0.68 


10.7 
15.2 


CQl 




5-10^ 














CQ2 


0.03 


10^ 


0.245 


-160 


66 


0.3 


0.71 


11.74 


CQ3 




102 






152 


3.6 


0.59 


11.17 


WRl 


0.03,0.04 
0.05,0.06 


50 


0.245 


-160 


203,294 
351,396 


1.1,1.8 
2.3,2.6 


0.68,0.54 
0.5,0.42 


11.35,10.60 
10.15,9.67 


WR2 


0.03,0.04 
0.05,0.06 


100 


0.245 


-160 


152,220 
266,302 


0.8,1.3 
1.8,2.0 


0.59,0.53 
0.49,0.39 


11.17,10.44 
9.80,9.33 


WR3 


0.03,0.04 
0.05,0.06 


200 


0.245 


-160 


102,150 
182,206 


0.5,0.9 
1.2,1.4 


0.57,0.50 
0.44,0.38 


11.10,10.40 
9.80,9.27 



Table 1: Parameters and characteristics of the dynamo models: the a-effect parameter Cq,; 
parameter controls the helicity dissipation rate; lo-^Qi^e ratio between the max- 

imum background turbulent diffusivity and the value at the top boundary (in the reference 
model of Stix (2002) this value is 0.245); tAJ^]^ is the minimum of the density gradient 
height at the top; (S^^)^^ is the maximum strength of the toroidal magnetic field averaged 
in the range of 0.95 — O.99i?0; B^^^^^ is the maximum strength of the radial polar magnetic 
field at O.99i?0; Cw is a calibration coefficient used for calculation of the sunspot number 
parameter; Aw is the sunspot number asymmetry parameter obtained in the models; Period 
is the period of sunspot cycles from the dynamo models. 



3. Results 



We summarize the parameters and characteristics of the dynamo models in Table 1. 
For simulating the sunspot number obtained from solar observations we use the relation 



suggested by Bracewell (1988): 



W{t) = Cw max(|(5(r,^,t))o95_o99^| ,0 < 9 < 180)'^' (3) 

where B is the toroidal magnetic field strength and Cw is the calibration coefficient, we 
use Cw = 1/40. We define the asymmetry parameter of the cycle as the ratio between the 



^^^^^ . The simulations 



dtW 

modulus of the mean decay rate and mean rise rate, Aw = 

\dtW\ 

' ' dtW<0 

were started from initial states with zero toroidal magnetic field and weak poloidal field 
which is symmetric about the equator. The solution is found by a semi-implicit method 
using a finite-difference approximation in radius and a pseudospectral decompisition in 
terms of Legendre polynomials in latitude. The numerical scheme conserves the parity 
of solution with respect to the equator. The characteristics of the dynamo models were 
determined from the stationary periodic solutions. 



3.1. Effects of the near surface diffusion and turbulent pumping 

In this part of the paper we fix the a-effect parameter Ca = 0.03 (the dynamo 
instability threshold is ^ 0.02 ). In the near surface layers the Coriolis number is rather 
small. Therefore, the turbulent pumping primarily depends on the density gradient A^^^ 
and diffusion coefficient tjt- The gradient parameter rA^^^ varies from ^ —7 at r = O.Tli?© 
to — 160 at r = 0.99/?©. To illustrate the inffuence of the pumping effect on the dynamo 
model solution we examine three different cases (PI, P2 and P3 in Table 1). Model PI 



employs the density gradient profile provided by the Stix (2002) model. In model P2, we 
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introduce an artificial limit on the level of rA^^^ = — 30 suggested by Kit chat inov et al. 



(2000). In model P3, we completely neglect the pumping effect. Profiles of the radial and 
the latitudinal pumping velocities at = 45° for models PI and P2 are shown in Figure Id. 
Note, that compared to the plotted values the amplitude of the velocities in the models is 
reduced by factor = 0.05. 

The time - latitude toroidal magnetic field "butterfly" diagrams, which were averaged 
over the depths from r = 0.95/?© to r = 0.99/?©, and the radial magnetic field evolution 
at r = O.99i?0 for models PI, P2 and P3 are shown in Figures 2. We find that the larger 
amplitude of the downward turbulent pumping results in the greater strength of the near 
surface toroidal magnetic field. The turbulent pumping increases the efficiency of the 
subsurface shear generation effect. This leads to a faster migration rate of the toroidal 



magnetic field to the equator (see, Yoshimura 1975). 



A similar effect can be produced by changing the turbulent diffusivity profile near the 
surface. In Figure 3 we show the results for two cases of the increased (Dl) and decreased 
(D2) turbulent diffusivity (Figure Ic). The diagrams for these two cases can be compared 
with the reference Stix's model PI in Figure 2a. The smaller the surface turbulent diffusivity 
level, the greater toroidal magnetic field strength. In the reference case the typical magnetic 
field strength is ^ 1 kG. If the surface turbulent diffusivity is three time smaller the field 
strength increases to ^ 1.6 kG; if the diffusivity is two time larger, then the field strength 
is only about ^ 450 G. The inclination of the toroidal field patterns migrating toward the 
equator (latitudinal migration speed) does not change considerably with the changes of the 
turbulent diffusivity profile. However, when the surface turbulent diffusivity is smaller, the 
toroidal magnetic field migrates closer to the equator, and the cycle becomes longer. 

For simulating the sunspot number obtained from solar observations we use the relation 
motivated by Bracewell (1988), see Eq.([3]). The estimated sunspot number for the models 
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with the increased and decreased sub-surface turbulent diffusivity is shown in Figure 3 
(left panels). The asymmetry between the rise and decay phases is clearly seen for the 
model with the decreased surface turbulent diffusivity. In the mean-field dynamo concept 
the decay phase of the large-scale toroidal magnetic field is defined by turbulent diffusion 



(Parker 1979). The decrease of the surface turbulent diffusivity increases the decay time 
when the toroidal field is located closer to the surface. 



3.2. Magnetic helicity effect and the Waldmeier's relations 

The evolution of magnetic helicity based on the conservation law is described by Eq.Q. 



Without helicity fiuxes from the dynamo domain (Kleeorin et al. 2000; Vishniac & Cho 



2001) or in absence of helicity dissipation, the generation of magnetic helicity by dynamo 
leads to "catastrophic quenching" of the a-effect, which stops the dynamo process (see, e.g., 
Vainshtein & Cattaneo 1992; Kleeorin et al. 2000; Brandenburg fc Subramanian||2005 ). In 



our model the dissipation of magnetic helicity is described by a decay term: — — — . We 
illustrate the catastrophic quenching in dynamo model CQl in Figure 4, which shows the 
results for i?^ = 5 x 10"^. In this case, the rate of the helicity loss from the Sun is small, and 
the dynamo process stops after 3-4 periods (2 magnetic cycles). The time-latitude diagram 



for the current helicity, which is estimated as = b-V x b x/^^? is shown together 
with the magnetic butterfiy diagram. We see that the total magnetic helicity generated 
by the dynamo is decaying much slower than the dynamo waves. Figure 4 also shows the 
evolution of the sunspot parameter, the total magnetic fiux, the total turbulent magnetic 
helicity and total large-scale magnetic helicity. For comparison with case i?^ = 5 x 10^ 
(model CQl), we show the results for = 10^, 10^ (models CQ2 and CQ3). In the case 
of = 10^ (model CQ2) the dynamo is stabilized at a quite low level with the maximum 
toroidal magnetic field strength of about 50 G and the total magnetic fiux of about 10^"^ 
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Mx. For the high dissipation rate of the magnetic hehcity, = 10^, (model CQ3) the 
maximum of the toroidal magnetic field strength inside the convection zone is about 450 G, 
and in the surface layer it reaches about 150 G. In model CQ3, the total magnetic fiux is 



about 4 X 10 Mx. This roughly agrees with observational results of Schrijver fc Harvey 



(1994). Taking into account the results shown in Figure 4d we can estimate the amount of 
the helicity loss from the Sun per cycle. In the model CQ3 it is about 1.26 x 10^^ Mx^. 

Model CQ3 is most relevant for comparison with observations. This case qualitatively 
reproduces the basic features of the solar cycle. Figure 5 shows snapshots of toroidal and 
poloidal fields, the time-latitude diagram of the near surface toroidal magnetic field, and 
the current helicity evolution, and the time-radius diagrams for the magnetic field and 
current helicity for this model. The time-laltitude diagrams illustrate the migration of the 
toroidal and poloidal fields and polarity reversal. The time-radius diagrams show migration 
of the magnetic field with radius at 30° latitude, and an interesting concentration of the 
field at r/i?0 ^ 0.9 — 0.92, or 60 — 70Mm below the surface. This concentration is related 
to the second maximum of the dynamo wave when it propogates from the bottom of the 
convection zone to the surface. 

Our results show that the current helicity changes the sign in the near surface layers 



at the begining of the cycle. A similar behaviour was found in observations of [Zhang et al. 
( |2010D (see their Fig.2). Our initial comparison of the current helicity pattern produced 



by the model reveals some disagreements with observations, e.g., the model does not show 
the change of the current helicity sign near the equator at the end of sunspot cycle. This 
problem needs a separate study, which should include a more sophisticated description of 
the helicity fiuxes. The time-radius diagram at latitude 30° for the toroidal magnetic field 
and the current helicity (Figure 5e) shows that that in the bulk of the convection zone the 
current helicity does not change much with time, and that the helicity is nearly constant 
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near the bottom of the convection zone. In the north hemisphere the magnetic hehcity is 
positive at the bottom of the convection zone because the kinetic part of the a effect is 



negative there (also see, Pipin & Kosovichev (2011) and Figure 1 there). In the near-surface 



layer the sign of the current helicity changes at the rising phase of the cycle. 

Model CQ3 clearly shows asymmetry between the growth and decay phases of the 
sunspot number parameter W{t). We find that the asymmetry increases with increase of 
the amplitude of the cycle, i.e., with increase of a-effect parameter Ca (see Figure 6). It 



is expected that for the higher Co, the dynamo period is shorter (see, e.g., Parkerj 1979) . 
This motivates us to look at the period - amplitude relationship for our model and also at 
the amplitude dependence on the growth/decay rate by computing a series of models WRl, 
WR2, WR3 for various values of Ca and R^. The results for the asymmetry parameter 
as well as the cycle period are summarized in Table 1. We compare the model results 
with the asymmetry estimated from the monthly smoothed sunspot number provided by 
the |SIDC[ The data set was additionally smoothed by means of the Wiener filter. After 



this, we divided the whole data set covering the time period from 1749 to 2010, into 
separate sunspot cycles. The cycles were divided by a program that catches the sequences 
of the sunspot minima. For each cycle we estimate the growth rate by a ratio of the cycle 
amplitude to the growth time. Similarly, the decay rate was defined. 

Figure 7 compares the model with these estimates in the form of the Waldmeier's 
(1935) relations: a) amplitude - rise rate, b) amplitude - decay rate, c) period - amplitude, 
d) rise time - amplitude, e) rise vs decay rates and g) rise vs decay times. The results 



obtained from the experimental data set confirm the findings of other authors (see, Vitinsky 



etaL][T986l [Hathaway et"aL][2002t [C ameron fc Schussler||2007[ [Kitiashvili fc Kosovichev||2009[ 



Karak & Choudhuri 2011). The computed dynamo models (WRl, WR2 and WR3) for a 



given range of the a-effect parameter = 0.03 — 0.06 and magnetic helicity dissipation 
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rate — 50 — 200) correspond reasonably well to the data points. However, there are 
differences. One possible source of the difference between the model and the data is clarified 
in Figure 7(e), which shows correlation between the rise and decay rates in the solar cycles. 
We find that for most solar cycles the rise rate is higher than the decay rate. The mean 
asymmetry parameter is Aw 0.68 ± 0.31 . As seen in Figure 7, our models have smaller 
A\Y . Thus, our dynamo models produce more asymmetric W profiles than the observed 
sunspot number. Increasing or decreasing the helicity dissipation by changing did not 
improve the agreement with the observations. 

It is clear that more studies of the turbulent properties of the Sun are necessarily. 
Nevertheless, the initial results of the dynamo models shaped by the subsurface shear layer 
are encouraging. 

4. Discussion and conclusion 

The paper presents a study of a solar dynamo model operating in the bulk of 
the convection zone with the toroidal magnetic field flux shaped into the time-latitude 
"butterfly" diagram in the subsurface rotational shear layer. We explore how this type of 
dynamo may depend on the radial variations of turbulent parameters and the differential 
rotation near the surface. The mean-field dynamo model takes into account the evolution 
of the magnetic helicity and describes its nonlinear feedback on the generation of the 
large-scale magnetic field by the o^-effect. We compare the magnetic cycle characteristics 
predicted by the model, including the cycle asymmetry and the duration - amplitude 
relation (Waldmeier's effects) with the observed sunspot cycle properties. We show that the 
model qualitatively reproduces the basic properties of the solar cycles. However, the model 
cycles are systematically more asymmetric than the observed cycles. 
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In section 3.1, it was shown that the radial profiles of the turbulent diffusivity and the 
density stratification scale in the sub-surface layer (between 0.95 — O.99i?0) may significantly 
infiuence the dynamo properties. In particular, the surface-shear shaped dynamo model 
favors the negative gradient of turbulent diffusivity in the sub-surface layer as follows from 
the standard solar model. This is contrary to the positive gradient of turbulent diffusivity 
often used in fiux-transport dynamo models (e.g., Karak fc Choudhurl||2011 ). In our model. 



a steeper gradient of the magnetic diffusivity results in a stronger toroidal magnetic field, 
a higher latitudinal migration speed and a longer magnetic field decay time in the surface 
layer. We found that the downward turbulent pumping of the horizontal magnetic field 
(associated with either toroidal or meridional magnetic field components) brings the dynamo 
properties in better agreement with observations, increasing the period of the magnetic 
cycle for a given turbulent diffusivity profile. The model shows the asymmetry between the 
rise and decay rates (and duration phases) of the toroidal magnetic field. The asymmetry 
increases with the increase of the turbulent diffusivity gradient in the sub-surface layer. 

The models shows a clear dependence of the asymmetry parameter (the ratio between 
the cycle's decay and rise rates) on the magnetic cycle strength. We compared a sunspot 



number parameter previously suggested by Bracewell (1988p with statistical properties 



of the solar cycle. Our model qualitatively reproduces the known properties, such as the 
Waldmeier's relations and the period - amplitude dependence. In particular. Figure 7e 
shows that the asymmetry is one of the basic features of the sunspot cycle activity (see 



also, Vitinsky et al. 1986). 



In our model this asymmetry depends on the parameters of the turbulent diffusivity 
in the near surface layer and on the rate of the magnetic helicity dissipation. If the 
magnetic helicity dissipation rate is higher, the asymmetry is smaller. The magnetic helicity 
dissipation rate infiuences the amount of the total magnetic fiux produced in the Sun. 



- 16 - 



According to Schrijver & Harvey (1994) the total magnetic flux produced during a solar 
cycle is about lO^^Mx. The models presented in the paper satisfy this constraint. The 
estimated amount of magnetic helicity loss in the dynamo model is about 10^^ Mx^ per 
cycle. 

Thus, we conclude that the dynamo models with the subsurface shear layer can satisfy 
the global constraints on the total magnetic flux produced by the dynamo and are able 
to qualitatively reproduce the known statistical properties of the solar cycle, like the 
Waldmeier's effects and the period - amplitude relation. We expect that the model can 
be further developed taking into account more accurately the turbulent properties of the 
sub-surface shear layer. The accurate description of the magnetic helicity dissipation is also 
important for the future progress. Direct numerical simulations of the convective turbulence 
and helioseismological data analysis techniques should help to improve our knowledge of 
the subsurface shear layer and the physics of solar dynamo. 
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6. Appendix 

We describe some parts of the mean-electromotive force. The basic formulation is given 

( H) 

in P08. For this paper we reformulate tensor Q^^j , which represents the hydro dynamical 
part of the a-effect, by using Eq.(23) from P08 in the following form, 



{^VT (e • A^'')) + (e • A^^)) ) } + (4) 
+ e.e,{3ry^(/r^(e.A(^))+/f (e-AW))} 
+ Svr { {e^A^^ + e,Ai^) ft^ + [e^Af + e, Af )) ft^ } . 

The contribution of magnetic helicity x = ^'b (a is a fluctuating vector magnetic fleld 
potential) to the a-effect is deflned as = C\'f^x-> where 



The turbulent pumping, 7jj, is also part of the mean electromotive force in Eq.(23)(P08). 



Here we rewrite it in a more traditional form (cf, e.g., Brandenburg fc Subramanian| 2005 



Radler et al. 2003) 



The effect of turbulent diffusivity, which is anisotropic due to the Coriolis force, is given by: 

r]ijk = 3r7r | (2/^''^ - /^"^^^ Sijk - 2f[''\ienenjk^ • (7) 

Functions f^^f\^^ depend on the Coriolis number f^* = 2tc^^ and the typical convective 
turnover time in the mixing-length approximation: Tc = i/u\ They can be found in P08. 
The turbulent diffusivity is parametrized in the form, tjt = Crjr]!^\ where r]!^"^ = ^ is 
the characteristic mixing-length turbulent diffusivity, is the RMS convective velocity, i 
is the mixing length, (7^ is a constant to control the intensity of turbulent mixing. The 



others quantities in Eqs.( 4[6[7 ) are: p is the density stratiflcation scale. 
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A(») = V log ( Ty^^") is the scale of turbulent diffusivity, e = Ct/ \ Q\ is a. unit vector along 
the axis of rotation. Equations (4][6|7) take into account the influence of the fluctuating 
small-scale magnetic flelds, which can be present in the background turbulence and stem 



Frisch et al. 


1975; 


Moffatt 


1978; 


Vainshtein 



& Kitchatinov 1983; Kleeorin et al. 1996; Brandenburg & Subramanian 2005|. In our 



paper, the parameter e = 



b2 



, which measures the ratio between the magnetic and 



kinetic energies of fluctuations in the background turbulence, is assumed equal to 1. This 
corresponds to the energy equipartition. The quenching function of the hydrodynamical 
part of a-effect is deflned by 



128;5^ 



(l6/3^ - 3 - 3 (4/5^ - 1) -^^^^ 



(8) 



Note, in notation of P08 i/ja = — 3/4(/)g^\ and /3 = — 



Lb 
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Fig. 1. — Parameters of the solar convection zone: a) contours of the constant angular 
velocity plotted in the range (0.75 - 1.05)^0 with a step of 0.025rio, ^^o = 2.86 • IQ-^s'^; b) 
turnover convection time rc, and the RMS convective velocity u'^] c) the background turbulent 
diffusivity 77^^ profiles; the solid curve shows the profile of the Stix(2002) model; the dashed 
and dotted curves show models with the reduced and increased sub-surface diffusivities 
discussed in Sec. 3.1; d) the radial (thick hues) and the latitudinal (thin lines) turbulent 
pumping velocities at 6^ = 45° for models PI (solid curves) and P2 (dashed curves), see Table 

1; 
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Fig. 2. — Illustration of the influence of the magnetic pumping effect (models PI, P2 and P3 
in Table 1). Time-latitude diagrams of the toroidal magnetic fleld averaged over the depth 
range 0.95 — 0.99Rq (gray scale) and the radial fleld at r = 0.99/?© (shown by contours), 
for three cases of turbulent pumping: a) model PI : the density gradient proflle provided 



by the Stix (2002|) model, the radial fleld varies in range ±21G; b) model P2: the density 
gradient effect is restricted as suggested by Kitchatinov et al. ( |2000[ ), the radial fleld varies 
in range ±14.5 G; c) model P3: the pumping effect due to the density gradient is neglected, 
the radial fleld varies in range ±1.4G. 
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Fig. 3. — Left column, the same as in Figure 2 for the larger (top) and smaller (bottom) sub- 
surface turbulent diffusivity, models Dl and D2 respectively. Right column - the simulated 
sunspot number, see definition in the text. 
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Fig. 4. — Illustration of the dynamo models with the catastrophic quenching (model CQl) 
and without it (models CQ2 and CQ3): a) toroidal field (gray scale) and current helicity 
(contours, in range ±3 10~^ G^m~^) for model CQl; evolution of global characteristics: b) 
the sunspot number parameter W, c) the total magnetic fiux, d) total large-scale magnetic 
helicity; e) total small-scale magnetic helicity for models CQl (dots), CQ2 (dashed curves), 
CQ3 (solid curves). 
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TIME, lYR] 



Fig. 5. — Illustration of dynamo model CQ3 : a) snapshots of the toroidal (gray scale) and 
poloidal (contours) magnetic field evolution for a half of the magnetic cycle; b) the time- 
latitude diagram for the toroidal magnetic field averaged over the depth range 0.95 — 0.99/?© 
(gray scale) and the poloidal field (contours, in range ±3 G); c) the time-latitude diagram for 
the toroidal magnetic field (gray scale) and the current helicity (contours, in range ±2 10~^ 
G^m~^); d) the time-radius diagram at latitude 30° for the toroidal magnetic field (gray 
scale) and the poloidal field (contours, in range ±3 G); e) the time-radius diagram at latitude 
30° for the toroidal magnetic field (gray scale) and the current helicity (contours, in range 
-2 10-4^2 10-^ G^m-i). 



""O 5 10 15 20 25 30 

TIME, [YR] 

Fig. 6. — The asymmetry of sunspot number parameter W in the series of dynamo models 
WR2 (Table 1) calculated for various values of Ca'- 0.03(dots), 0.04 (dot-dashed curve), 0.05 
(dashed curve) and 0.06 (solid curve). 
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Fig. 7. — The Waldmeier's (1935) relations in the models WRl (dashed line), WR2 (solid 
line) and WR3 (dot-dashed line): a) amplitude - rise rate, b) amplitude - decay rate, c) 
period - amplitude, d) rise time - amplitude, e) rise - decay rates, g) rise - decay times. The 
results from the SIDC sunspot data are shown by circles. 



